One-loop divergences in the Galileon model 

Tiberio de Paula Netto§ and Ilya L. Shapiro § 

Departamento de Fi'sica, ICE, Universidade Federal de Juiz de Fora 
Juiz de Fora, 36036-330, MG, Brazil 



Abstract. The investigation of UV divergences is a relevant step in better 
understanding of a new theory. In this work the one-loop divergences in the 
free field sector are obtained for the popular Galileons model. The calcula- 
tions are performed by the generalized Schwinger-DeWitt technique and also 
by means of Feynman diagrams. The first method can be directly generalized 
to curved space, but here we deal only with the flat-space limit. We show 
that the UV completion of the theory includes the 7rD 4 7r term. According 
to our previous analysis in the case of quantum gravity, this means that the 
theory can be modified to become superrenormalizable, but then its physical 
spectrum includes two massive ghosts and one massive scalar with positive 
kinetic energy. The effective approach in this theory can be perfectly suc- 
cessful, exactly as in the higher derivative quantum gravity, and in this case 
the non-renormalization theorem for Galileons remains valid in the low-energy 
region. 

MSC: 81T15, 81T18, 83D05 
PACS: ll.10.Gh, 04.50.Kd 

1 Introduction 

The Galileons are qualitatively new models of scalar field with numerous applications. 
Originally Galileons were introduced in the context of Dvali-Gabadadze-Porrati model of 
gravity [T] and attracted a great deal of interest (see, e.g., [21 [31 H] and further references 
therein). One of the standard motivations to consider Galileons is that they are assumed 
to possess some unusual renormalization properties [HI [6], especially when treated in the 
effective quantum field theory framework [TJ El [9] . 
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In short, Galileons are second derivative theories with much more derivatives in the 
action. This means that the kinetic term of such a model has only two derivatives, while 
the interacting terms have much more derivatives. As a result, the tree-level propagator 
of the theory is free from higher derivative ghosts and, at the same time, the possible 
divergences have so many derivatives that the terms which are present in the classical 
action, never get renormalized j5j E]. In this respect, the Galileon model is interesting 
to compare with the fourth derivative quantum gravity (HDQG) theory [TO] (see also 
[TT] for detailed introduction) which was further generalized in [12]. The common point 
between the two theories is the presence of higher derivatives. In case of HDQG this 
makes the theory renormalizable [10] or even superrenormalizable [12]. However, there is 
a price to pay: the particle spectrum of the theory includes higher derivative ghost pU] 
or ghosts [12]. At the same time, Galileon model is very much nonrenormalizable, in a 
sense that there are many possible divergences, however all of the possible logarithmically 
divergent terms have much more derivatives that the initial classical action, as a result 
the low-energy sector of the theory is free from strong UV quantum corrections. 

The described scheme of constructing a theory which is not affected, in the UV limit, 
by quantum corrections, looks very attractive, but there is one important point to verify. If 
there are higher derivative divergences in the propagator sector, then the UV completion 
of the theory actually has massive ghost. Eliminating this ghost from the spectrum 
leads to the possible unitarity breaking, exactly as in the HDQG. And, exactly as in the 
HDQG, we can try to formulate the theory in such a way that the ghost is not generated 
at relatively low energies. So, the program of exploring Galileons at quantum level should 
be supplemented by direct calculation of the quantum contributions to the scalar field 
propagator. Such a calculation, at the one-loop level, is the subject of the present paper. 

In general, the status of the quantum theory essentially depends on the structure 
of UV divergences in the free field sector. It may happen that the quantum theory 
is non-renormalizable, or it can be completed to become renormalizable or even super- 
renormalizable, like the HDQG model of |12j . In the last case the one-loop contribution 
may be the only one which is relevant. In what follows we shall figure out what is the 
situation, namely whether the Galileons theory is non-renormalizable, renormalizable or 
superrenormalizable and whether it can be unitary in a strong sense or only as an effective 
field theory. 

The paper is organized as follows. In Sect. |2] we present a very brief description of 
the classical action of the model. In Sect. [3] the derivation of one-loop divergences is 
performed by means on the generalized Schwinger-DeWitt technique [13] with expansions 
originally developed in [14J. The main advantage of this technique is that it can be also 
applied to the similar derivation in curved space, however in the present Letter we restrict 
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our attention only by the flat-space limit, which is indeed sufficient to answer the questions 
formulated above. In Sect. HI which was written for better control and for illustrative 
purposes, we show how the same result can be obtained by means of a more conventional 
Feynman diagrams technique. In Sect. [5] we present some additional discussions of the 
result and draw our conclusions. 

2 Brief description of the model 

In a four dimensional space-time there are only five Lagrangians with single scalar field 
7T, which are invariant (up to a total derivative) under the following transformation 

7T 7T + + C, (1) 

where c and are constants. The transformation (CEJ) is called Galilean transformation 
and the field tt is called Galileon. These five Lagrangians can be represented by the 
following structures: 

£1 = 7T, 

U = ^7r^7r(D7r) 2 - d^d^d^U-K - ^-d^ixd^d" nd e ixd p ix 

£ 5 = ^7r^7r(D7r) 3 - ^Trd^Tn^TrDTr - -d^Trd^d^d^d^nTt 

+ ~d fl d l/ Trd* 1 d u Trd e Trd e d X Trd x <K - d^d^d v Txd u d Q -nd Q d x d X Ti . (2) 

The full Lagrangian for the field n is a linear combination of the above Lagrangians 

5 

H-k = ^ ^ CiCi , (3) 

where Cj's are generic coefficients. 

3 Calculation of the one-loop counterterms 

In this section we shall present the details of the calculation of the one-loop counterterms 
of the theory (J3]). For the purpose of calculating the divergences we shall apply the 
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background field method (see, e.g., Chapter 2 of [TT] for introduction) and the generalized 
Schwinger-DeWitt technique [13]. Let us start with the usual splitting of the fields into 
background and quantum part 

7T — ► 7r' = 7T + a. (4) 

The one-loop effective action is given by the expression 

= ^Trlntf, (5) 

where H is the bilinear form of the action given by Lagrangian fl3]). Substituting flj]) in 
(J3]) one can find the bilinear form of the action 

= ~J d^x a Ha, (6) 



where H has a form 



H = □ + A + A + ■ ■ ■ , (7) 



here A ~ ^K 71 ")) A ~ CI 71 " 2 ) an d so on - Let us note that the scheme of calculation which 
we are going to perform is designed to get the correction to the propagator, therefore we 
do not need to take into account the terms beyond A and the square of the A-term in 
this and further expansions. 

In order to calculate the divergent part of the one- loop effective action (jSJ), in the 
second order in tt, one can perform the expansion, which is similar to the one which was 
previously used in p3j (described also in details in Chapter 8 of the book [TT]). 

Tr InH = Tr In (□ + A + A H ) 

= Tr In □ + Tr In f 1 + A— + A- + • • • ) 

Van / 

= Tr lnD+ Tr (A- + A- - -A-A-) + ■•■ (8) 

The omitted terms are 0(tc 3 ) and hence (as we have already mentioned above) they are 
actually irrelevant for our purposes. In order to reduce the amount of calculations we 
shall consider only the particular case of the flat background space-time. Then the only 
type of universal trace that does not vanish has the form 



Tr dn, ■ ■ ■ d ur> . — 



2i r a (n ~ 2) 



div ~ e J 2«-2(n-l)!' 1 ' 
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with n > 2 and 



(n-2) 



W2 #m 3W • • • ^2n- 3 ^2„-4 + all permutations . 



(10) 



By using the formula in we can see that the trace of the p2(l/n)-term corre- 
sponds to n — 1 in Eq. fl9]) and hence it is finite in dimensional regularization. Therefore, 
only the last term of equation (jSJ) gives contribution to the divergences. Then, 



Tr lnPT 



div 



-ilr A-A- 

2 □ □ 



div 



To calculate (II ip we need an explicit form of Pi, namely 



Pi 



d^d v , where 



4c, 



(12) 



The commutations which enable one to reduce the problem to the universal trace (Q 
are performed as follows: 



Tr A — A — 

□ □ 



div 



div 



Tr U^d^lu^d.d^ + 



Furthermore, the last commutator can be transformed as 



□ 



d, L d u 



1 

□ 



(/ill 



(13) 



+4(«9 A nf>HcW^ " ^(d^W^A-W^ - 2(9 A f>^)9 A ^^ 



+ 16 (^d p d x d T U^) d w d p d x d T d v d^ + [je)- 

The terms with background dimension (we assume the reader is familiar with the termi- 
nology of p2]) of more than can be safely omitted here because they do not contribute 
to divergences. Substituting the equation ( IT4l) into ( TT31) we finally get, after some tedious 
calculations, the following expression: 



Tr A-A- 

□ □ 



div 



+ 8U^(d a d x nu^d x d,d^ - i2U^{d x d T nu^)d a d p d x d T d,d u ^ 
+W aP {d a d p d x d T U^)d x d T d^d u ^ i - lW aP {d a d p d x d T U^)d p d p d x d T d v d l 1 



M D 5 



+ lW<*P(d fJ d p d x d T U>* u )d a dpd 0J d p d x d T d l ,d fl — 



div 



(14) 



Now we are in position to use the universal trace (Q. In this way, after certain work, 
we obtain the following traces [here e = (47r) 2 (n — 4) in dimensional regularization] : 

1 



- Tr U^idadpDU^d^ 



□ 3 



div 2,6 



d\ U a %d a dpnufi, 



Tr U ap {u 2 u^)d a dpd„d v 



□ 4 



div 



12 e 



8Tr U a '\d a d x UUndpd x d,d v —^ 



div 



2i 
3e 



d 4 x \ 2U^(d a d^nu fM/ ) + u a/3 (d a d p nu£) 



4Tr U al5 (d a dpd x d T U^)d x d T d,d v ^ 



div 



3e 



d 4 x \2 U a ^{d a d p d ft d v U fa ') + U a)3 {d a dpnun 



- 12 Tr U a ^(d x d T aU^)d a d^d x d T d fl d u 



div 



= ^ y d 4 x {u^(d a d^aup + ~ f>:(n 2 f/;) + 4U:(d a d,nu»n 



and 



16 Tr U^{d a d p d x d T U^)dpd p d x d T d v d, — 



16 Tr U aP (d J} d p d x d T U» v )d a dpdud p d x d T d v d lx \ 



div 



+ 2U a \d a d^d v iin + 4£C (c^n^) + ~ fV (n 2 ^) } . 

By using the traces listed above, Eq. (fl4"|) can be reduced to 



Tr A - A - 

□ □ 



(/in 



| y rf 4 x { 1 u^idjpnup - i ^(n 2 ^) 



(15) 



(16) 



Replacing, term by term, the explicit form of operator given by Eq. (|T2|) into Eq. 
(|T6l) . after some algebra we arrive at 

2i 



Tr A — A — 
□ □ 



-d \ (fx 7rn 4 vr. (17) 
div e " ' 



Finally, from this expression and the equations flSJ), ( ITTj) we obtain the result for the 
one-loop divergences of effective action in the free field sector, 



r£ = -f| y^TrnV (is) 

The expression (|T8|) is providing us some relevant information about the quantum 
properties of the Galileons theory. Let us present it in the systematic form. 

• The UV logarithmic divergences of the theory require 7rD 4 7r-type counterterms. This 
means that the consistency of the theory requires that the same term is included into the 
classical action of the theory. If we do not include such a term, it will emerge with infinite 
coefficient anyway and then no control over this term via the effective approach will be 
possible. 



• According to the analysis of the similar gravitational theory in [12], the theory with 
an extra 7rD 4 7r-type classical term has, typically, two massive ghosts (excitations with 
negative kinetic energy), one massive scalar degree of freedom and, of course, the "orig- 
inal" massless scalar mode. This situation means that the propagator of the field 7r has 
the general structure 



I A 1 , A 2 A ± 

- Ti 

^1,2,3 > , m 3 > m 2 > mi > . 



G{P) = -n 5^ 2 + ~2~, 2 2^ 2' ( 19 ) 



Here we assumed Euclidean signature and that there are no tachyons in the spectrum. 
The last can be always provided by adjusting the coefficients of subleading 7rD 2 7r-type 
and 7rn 3 7r-type terms. Let us note that similar, curvature-dependent, terms are likely 
to be requested also by the divergences in curved space-time. On the other hand, all 
considerations presented below should be valid also in the presence of tachyons. 

• According to the analysis of the similar gravitational theory in [12J, the theory with an 
extra 7rD 4 7r-type classical term is superrenormalizable, for any choice of the coefficients Ck 
in the interacting sector of (T5]). This means, in our case, that the divergences can emerge 
only at the one-loop level. Starting from the second loop only one-loop sub-diagrams 
can be divergent. This feature does not depend on the flatness of space-time and will 
definitely hold also in curved space-time case (see [TT] for introduction to the general 
theory of renormalization in curved space-time). 
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• The effective approach in the quantum Galileon theory is perfectly possible if we put 
a sufficiently small parameter in front of the classical J 7rD 4 7r-term in the action. Such 
parameter should have the form of M -6 and, therefore, the choice of a small coefficient 
of the higher derivative term means we choose a huge mass parameter M. The masses of 
both massive ghosts m 1; m$ and of the massive scalar with positive kinetic energy, 7712, 
will be of the same order of magnitude as M (see [12] for details). If we consider classical 
or quantum phenomena at the energies much smaller than M, the massive modes of the 
scalar do not become active and the conclusions of [3 [6] and (7J[8] remain correct, including 
the non-renormalization theorem. In the effective framework the "correct" quantum result 
is supposed to be the one we have derived above, and not the one of the complete theory 
with the J 7rn 4 7r-term. Finally, this means that ([T]) is the low-energy symmetry, so the 
fact that it is violated by the one-loop divergence (fT8l) is irrelevant in this framework. 

4 One-loop divergences from Feynman diagrams 

In the previous Section we calculated the one-loop logarithmic divergences by means of 
the generalized Schwinger-DeWitt technique. This approach has an advantage because 
it is relatively easy to generalize the results for a curved space-time. However, since in 
practise we are dealing only with flat space-time limit, here we perform an equivalent 
calculations by a more traditional Feynman diagrams-based calculation. The purpose of 
this section is to have an extra control of the result and also for the illustrative reasons. 
Since the theory looks unusual, this consideration may be instructive. Without going 
into full details, we will also provide a comparison between diagrams and universal traces, 
considered in the previous Section. 

Let us considerer the diagrams which contribute to the two-point function of the 
Galileon field at the one-loop level. The first set of diagrams which are all generated by 
the Lagrangian £3 is shown in Fig. 1. 




Figure 1. The first set of diagrams coming from the /Vterm vertex, 
which contribute to the propagator of the Galileon field. Here primes indicate 
derivatives. 

We will be interested in the behavior in the large momentum regime. The integral 
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associated with the first diagram from Fig. 1 is 

nib) = c 



2 / d A q p A q ll g u {q f " - p^){q v - p v ) 

3 ' (2tt) 4 q\q-V? 



(20) 



Expanding the denominator of equation f[2"Uj) . we found at large momentum scale the 
relation (in this equation we omit all tensorial indices and coefficients, and are only 
interested in power counting related to the divergences and powers of momentum p) 



ni(p) cl I dql ■■■+p 6 q + p 7 + — + 

Jo v q 



(21) 



One can easily see that this diagram contains ultraviolet logarithmic divergences and the 
power of p corresponding to this divergence is eight, which fits the D 4 -term obtained in 
the previous Section. 

The others diagrams of our interest are the ones shown in Fig. 2 




— 1 — y- 




</ 

v — J — y- 





V ' 7^ 



Figure 2. The second set of (tadpole-like) diagram provided by £3, which 
contribute to the two-point function. 



and also the graphs shown in Fig. 3, 




Figure 3. The diagrams generated by lagrangian £ 4 which may con- 
tribute for propagator of Galileon Model. 

In fact, the diagrams from Figures 2 and 3 do not contribute to the divergences in 
the case of Galileon theory. The reason is that these diagrams include derivatives of the 
propagator in a single space-time point, and hence vanish. 

Using this simple analysis of the diagrams one can explain qualitatively the contri- 
butions for the divergences for each term in the expansion for effective action ()8]). First 
let us consider the term TrPx^. Since P\ is 0(tt 3 ), it contains only contributions of 
the Lagrangian £3 and is proportional to C3. Then, TrPi^ is also proportional to C3 and 
we can see that this term does not contribute to divergences because it corresponds to 
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Feynman diagrams which are proportional to C3 and are zero after we use Wick's theorem, 
since £3 has an odd number of fields. 



then the contribution of this term are given by tadpoles diagrams of Fig. 3. As we have 
mentioned above, this diagrams do not contribute, hence we can see why this term makes 
no contribution to the counterterms. Finally, the last trace is TrPi^Pi^j. This term is 
proportional to c|, its contribution to the logarithmic divergences is different from zero 
and is given by diagrams of Fig. 1 . 

The considerations presented above enable us to write the expression for the divergent 
part of the two-point function as presented in Fig. 4. 



GSr rt (*, V) = i (a, x r+Ct+i + a, x ; + «3 x 



Figure 4. Diagrammatic representation of the Green function. 01,2,3,4 
are combinatorial coefficients. 

The definition of the full polarization operator is given in Fig. 5. 



The next term is TrP 2 ^. Remember that P 2 is 0(ir 4 ) and is proportional to c 4 , 





Figure 5. Full polarization operator. 



with 




(22) 




(23) 




(24) 



and 




(25) 
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To evaluate these integrals we used dimensional regularization, reformulating the the- 
ory in the space-time of 2cj complex dimensions, where the integrals are convergent. Using 
the formulas of the Appendix, we found the result 

nifow) = ^p 8 /i, 

Ui(p,u) = 0, 1 = 2,3,4. (26) 

The integral I\ is defined in the Appendix. To find the divergent part of the polarization 
operator we consider the limit u — > 2. The divergent part is given by the pole of r(2 — u) 
in I\. We find 

J iL = -|. ( 27 ) 

where e = (4ir) 2 (n — 4). Taking into account the combinatorial coefficient ai = 4, we 
arrive at 



IWp,u;) = '-fp'ir = - l -f P 8 . (28) 



The last part is to calculate the diveregent part of Effective Action in the coordinate 
representation. For this end one can use the Dyson's formula 

G- 1 = G, 1 - n(p) . (29) 

Then the divergent part of Effective Action can be written as 

= -m div (n) = -^n 4 . (30) 

07C071 e 

The equation f l30]) can be easily solved and we arrive at 

r£ = - c f e /^ttdV, (si) 

that is exactly the result obtained in previous section. 



5 Concluding discussions 

We have developed the background field method and calculated the one-loop divergences 
for the Galileon model. It turns out that the UV completion of the theory includes higher 
derivative sectors, as it was indeed anticipated in [51 |6]. An interesting new aspect is 
that this UV completion leads to the superrenormalizable quantum theory, where only 
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the one-loop contribution to the effective action is divergent and everything beyond the 
one-loop order is finite. 

There is an interesting similarity between the quantum Galileons model with this 
higher derivative completion and the higher derivative quantum gravity (HDQG). In fact, 
the unique conceptual difference is that the Galileon model with an extra J 7rD 4 7r-term 
is strongly superrenormalizable, while HDQG admits different levels of renormalizable 
and superrenormalizable theories. At the same time, the status of ghosts in these two 
theories is very close. In case of HDQG the Planck mass Mp plays the role of the massive 
parameter M which was discussed at the end of the section |3j In both cases one can 
provide the absence of ghosts at the tree level for sufficiently low energies. 

One simple test of the last statement has been applied recently in the cosmological 
framework [18] . It was shown that the physically relevant cosmological solutions in the 
higher derivative gravity theory (even with complicated semiclassical corrections) are 
stable with respect to graviton perturbations (gravitational waves). Definitely, this output 
is expected to hold only until we do not start to deal with the perturbations with the 
initial amplitude of the Planck order of magnitude. However, after the Universe passed 
through its initial Planck-scale epoch, such violent perturbations are never generated, 
and therefore the theory is safe at the classical level. Probably, this should mean that the 
quantum theory is also free of the ghost problem at the tree level. Of course, this is not an 
obvious statement, because it is not really clear how the most relevant classical solutions 
(such as cosmological one, for example) of the gravitational theory can be reproduced 
via the linearized gravity approach. At the same time, the stability of the cosmological 
solution [18J is definitely more fundamental issue than our skills in linearizing gravity, so 
we can definitely say that we have a strong positive arguments in favor of higher derivative 
theories^]. 

Let us, finally, discuss some practical lessons which we can learn from the analogy with 
the HDQG case. As far as the main applications of Galileons is related to cosmology, 
it would be definitely interesting to consider the stability of the classical cosmological 
solutions in the presence of higher derivative terms which are a necessary UV completion 
of the theory. For this end one has to complete the J 7rn 4 7r-term derived here by the 
corresponding curvature-dependent terms. Regardless of most of the relevant information 
for such a cosmological application can be perfectly well obtained from power counting 
arguments, it would be anyway reasonable to generalize our calculation of the one-loop 
divergences to the curved space case. 

3 There are many other interesting proposals towards the solution of the ghost issue in HDQG [THl 
[20l [2T1 l22l [23] which can be also productive. In any case it is important to care about higher derivative 
terms in gravity, since they are requested by consistency of the quantum theory of matter fields 124) . 
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6 An extra note on more general scalar-tensor theo- 
ries 



It is possible, in principle, to generalize the results considered above to the more compli- 
cated and general scalar-tensor model with second-order field equations. Such a theory 
has been recently considered in [IT] (see further references therein). The action of this 
general model does not satisfy the symmetry ([1]) and this opens the way for infinitely 
many new terms in the Lagrangian. The form of the curved-space Lagrangian is 

5 

C = ^2&, (32) 

i=2 

C 2 = K((f>,X), where X = -~c 
C 3 = -G 3 (0,X)D0, 

U = G^,X)R+ G 4iX [(D0) 2 -(V M V^)(V^ 
£ 5 = G 5 (<j ) ,X)G tlv - ^G 5 ,x[(n0) 2 -3(D0)(V / ,V^)(V^V ! ' 
+2(p<f>) (V»V a <p) ( V a V^) ( V" V, 



where Gk{4>,X), with k = 3,4,5, are arbitrary functions and Gkx{<f>,X) are the corre- 
sponding derivatives with respect to X. 

The question in which we are interested in is whether and how the result (ITS]) gets 
modified in the more general model (132]) . In order to address this issue, it is sufficient to 
perform the analysis of the power counting. Hence, our consideration of this model will 
be brief, so that we leave the details of the performed analysis as exercise for an interested 
reader and give only the main result. 

As far as we are interested only in the counterterms which contribute to the propa- 
gator of (f) and do not intend to quantize metric, the modifications which come from the 
curvature-dependent terms in (13"2"]) are irrelevant. Furthermore, according to our previous 
analysis, only the contributions to the vertices with three legs are significant, while the 
vertices with four and more legs play no role. Taking these two observations into account, 
one can easily see that the result strongly depends on the presence of a constant term in 
the function G^x- 

Let us assume that G 5 (<p,X) can be expanded as 

G 5 (0, X) = G 5O (0) + GM)X + G 52 (0)X 2 + ... , 
where G 51 (<f>)X = G 5 io + G 511 + G 512 2 + G 513 3 + .... (33) 

It is easy to see that if the coefficient G510 in ( 133]) is zero, then the result ( |T8]) does not 
change (except the coefficient, of course). However, in case G510 7^ the consideration 
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of the superficial degree of divergences of the relevant diagram shows that the leading 
counterterm will be different from (118jl . In this case we can expect the counterterms of 
the form 



where the dots indicate the presence of possible terms with lower powers of □. The higher 
order in derivatives in ( |34l) compared to ( fi~8l is because the Gsio-term in ( |32l) leads to 
the vertices with three legs and five derivatives, which are not present in the Galileon 
case ([2]). Qualitatively, the consequences of the higher derivative terms remain the same 
as it was discussed above. The theory with the corresponding UV completion would be 
superrenormalizable, and possesses a (larger, in this case) set of ghosts and massive scalar 
states with positive kinetic energy. 

Appendix. Massless integrals 

To calculate the integrals of Feynman diagrams from Sect. 4, we need the divergent parts 
of some massless integrals, which are given below. The basic formulas (!35|) - (l39|) can be 
found in [15] and other integrals can be obtained by the method explained in [16J. All 
integrals are defined over Euclidian space and must be understood though the prescription 
for massless case explained in [15] 



div 



(1) 




(34) 




(35) 




(36) 




(37) 




(38) 



/ 



d 2u3 qq^quq a qi3 



(39) 



(27r) 2 <V(g - p) 2 



d 2w q q^quq a qpq P 
(2TT)^q 2 (q-p) 2 



PfiPvPaPfiPpho + K^afSpIu + Lnvapplu , 



(40) 
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(2n) 2ul q 2 (q — p) 2 ~ P^PvPaPfiPpPu ' n n^al3pu: J 14 ~r dfivappu J 15 
"I - 'Fpuaffpui-^IG i 



where 



-E'/ii/a = ^i/Pa + all permutations , 

= S^papp + + all permutations , 

HfivaP = + a ^ permutations , 

K^app = bpvPaPfiPp + all permutations , 

£^ Q( 3p = 5pv5 a pP P + all permutations , 

RpvaPpu, = SpvPaPpPpPu, + all permutations , 

Stwappu, = SpvbapPpPu, + all permutations , 

Tpvotppu = S/xvSapSpoj + all permutations . 
The integrals J 2 , ... , he can be expressed in terms of the basic integral I 1: 

h = — 7 r r(2 - u)F(u - l)F(u - l)p 2 ^ 

(47r) w r(2w - 2) v ' K ' K ,F 



and are given by the expressions 



h = \h, 



T - ~ p2 J 
J 3 — 7777 TT-'l ) 



h = 77777 77 J l > 



4(2w- 


-1) 






2(2cu- 


-1) 


(w + 


1) 


4(2w- 


-1) 



j — 
15 



= 8(2^1) (55) 



_ (h, + l)(o, + 2) 
/7 " 4(4c 2 -l) Jl ' (56) 



-(u; + l)p 2 

/s = 8(4c 2 - 1) Jl ' ^ 



J 9 = TT J i . (58) 



16(4w 2 - 1)' 



(^+3)^+2) 

Jl ° = 8(40,2-1) Jl ' (59) 



-(a, + 2 y 

16(4w 2 - 1) V ; 



4 

/l2= 32(4^-l) 71 ' (61) 



_ (u; + 4)(a; + 3)(a; + 2) 
/l3 " 8(2c + 3)(4^-l) Jl ' (62) 



_ -(a, + 3)(a, + 2) P 2 
/l4 " 16(2c + 3)(4c 2 -l) Jl ' (63) 



'15 



'16 



32(2w + 3)(4o; 2 - 1) 



A , (64) 



p 6 

/i • (65) 



64(2w + 3)(4cu 2 - 1) 
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